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We report the first observation of B — > ip(2S)ir , a Cabibbo- and color-suppressed decay. This 
analysis is based on 657 x fO 6 BB events collected at the T(4S) resonance with the Belle detector 
at the KEKB energy-asymmetric e + e~ collider. The measured branching fraction is (2.44 ± 0.22 ± 
0.20) xl0~ 5 and the charge asymmetry is A = 0.022 ± 0.085 ± 0.0f6. The ratio of the branching 
fractions B(B ± -> ip(2S)-K ± )/B(B ± -> V(25')A' ± ) = (3.99 ± 0.36 ± 0.17)% is also determined. 



PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Gx 



The study of exclusive B meson decays to charmo- 
nium has played an important role in exploring CP- 
violation [lj. Among these, Cabibbo-suppressed charmo- 
nium decays provide an opportunity to understand non- 
leptonic B decays @, S 0, II @, 0, B Q and can be studied 
at the B factories which, due to their high integrated lu- 
minosities, overcome the suppression factor. At quark 
level, the decay B~ — > tp(2S)TT~ [lfj proceeds primar- 
ily via a b — > ccd transition; the leading contribution to 
this decay comes from the Cabibbo- and color-suppressed 
tree diagram shown in Fig. 1. The measurement of the 
branching fraction for the decay B~ — > ip(2S)n~ , which 
has not been observed so far, is important for detailed 
studies of the b — > ccd transition. 

Assuming tree dominance and factorization, the 
branching fraction B(B~ — > ip(2S)Tr~) is expected to be 
about 5% of that of the Cabibbo-favoured mode B~ — > 
tjj{2S)K~ [ll| . Under this assumption, the charge asym- 
metries for the B~ -> ip(2S)w~ and B~ -> ip(2S)K~ de- 
cays are expected to be negligibly small in the standard 
model. However, the penguin amplitude in the b — > ccd 
transition contains a complex phase due to one of the 
Kobayashi-Maskawa matrix elements [12j, Vtd, and there- 
fore, if penguin or new physics contributions are substan- 
tial, direct CP violation may occur in B~ — > yb(2S)ii~ 

In this paper, we report the first observation of the 
B~ — > i/j(2S)tt~ decay, along with a measurement of the 
branching fractions B(B~ — + i/j(2S)tt~) and the ratio 
B(B~ -> tjj(2S)tt-)/B(B- -> ip(2S)K~). A search for 
direct CP- violation in B~ — > ip(2S)ir~ decays is also pre- 
sented. These measurements are based on a data sample 
that contains 657 x 10 6 BB events, collected with the 
Belle detector [15j at the KEKB [16( energy-asymmetric 




B 



FIG. 1: Leading-order tree level diagram for the decays under 
study. 



e + e~ collider operating at the T(45') resonance. 

The Belle detector is a large solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD) 
surrounded by a 50-layer central drift chamber (CDC), an 
array of aerogel threshold Cherenkov counters (ACC), a 
barrel-like arrangement of time-of-flight (TOF) scintilla- 
tion counters, and an electromagnetic calorimeter (ECL) 
comprised of CsI(Tl) crystals. All these subdetectors are 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return yoke 
located outside the coil is instrumented to detect K a L 
mesons and to identify the muons (KLM) . The data sam- 
ple used in this analysis is collected with two different de- 
tector configurations. The first sample of 152 x 10 6 BB 
events is collected with a 2.0 cm radius beam-pipe and a 
3-layer SVD, while the remaining 505 x 10 6 BB events 
are collected with a 1.5 cm radius beam-pipe, a 4-layer 
SVD and a smaller-cell inner drift chamber 17J . The de- 
tector is described in detail elsewhere [15j|. A GEANT- 



bascd Monte Carlo (MC) simulation is used to model the 
response of the detector and determine the efficiency of 
the signal reconstruction [18|, |l9( ■ 

The reconstruction of the ip(2S) meson is performed 



using the £ + £~ (£ = e or fj) and J/tpir + ir~ decay chan- 
nels. The J/ip mesons are reconstructed in the £ + £~ 
decay channel. Both daughter tracks of J/ip — ► £ + £~ or 
ip(23) — > £ + £~ decays are required to be positively identi- 
fied as leptons. From the selected charged tracks, e + and 
e~ candidates are identified by combining information 
from the CDC (dE/dx), E/p (E is the energy deposited 
in the ECL and p is the momentum measured by the SVD 
and the CDC), and shower shape in the ECL. In addi- 
tion, ACC information and position matching between 
electron track candidates and ECL clusters are used in 
the identification of electron candidates. Identification of 
(i candidates is based on the track penetration depth and 
hit pattern in the KLM system. 

In J/iJj — ► e + e~ and ip(2S) — ► e + e~ decays, the four- 
momenta of all photons within 50 mrad of each of the 
original e + or e~ tracks are included in the invariant 
mass calculation [hereafter denoted as M e + e -( 7 )], in or- 
der to reduce the radiative tail. The reconstructed in- 
variant mass of the J/ip candidates is required to satisfy 
-0.150 GeV/c 2 < M e+e - (7) 
or -0.060 GeV/c 2 < M M+M - 
where rnj/^ denotes the nominal world-average J/ip mass 
[2fj|. These intervals are asymmetric in order to in- 
clude part of the radiative tails. In ip(2S) reconstruction, 
the invariant mass is required to satisfy 3.55 GeV/c 2 < 

,.- < 



m J/v , < 0.036 GeV/c 2 
mj/^, < 0.036 GeV/c 2 , 



AL 



'(7) 



< 3.75 GeV/c 2 or 3.65 GeV/c 2 < M 



3.75 GeV/c 2 . For the ip(2S) -> J/ipn + n~ candidates, 
AM = M i + e - 7T + 7T - — Mi+i- should satisfy the condition 
0.578 GeV/c 2 < AM < 0.598 GeV/c 2 . In order to re- 
duce the combinatorial background from low-momentum 
pions, M n + V - (invariant mass of the two pions coming 
from the ip(2S) decay) is required to be greater than 
0.40 GeV/c 2 . In order to veto peaking background com- 
ing from K% (B~ — » J/ipK*~(K^Tr~)), all pairs of pions 
with mass 0.4856 GeV/c 2 < M„+„- < 0.5096 GeV/c 2 are 
rejected in the ip(2S) — > J/ipir + ir~ reconstruction of the 
B~ — > ip(2S)n~ mode. A mass- and vertex-constrained 
fit is performed to all the selected ip(2S) and J/ip candi- 
dates in order to improve the momentum resolution. 

We combine ip(2S) and ir~ (K~) mesons to form B 
meson candidates. The combined information from the 
CDC (dE/dx), TOF and ACC is used to identify tt" 
(K~) on the basis of the ir — K likelihood ratio, TZ^^k) = 
£ir(K) / (£■* + Cr), where £„ (Ck) is the likelihood of a 
pion (kaon) hypothesis. Charged tracks with 7Z, n > 0.85 
(TZk > 0.85) are identified as 7r~ (K~). This requirement 
is 88.0% (82.6%) efficient for n (K) with a K (tt) fake 
rate of 6.0% (3.7%). 

To discriminate the signal from background, we use 
two kinematic variables: the beam-constrained mass, 
A4c = x/^bcam _ P*b &n ^ tne ener gy difference, AE = 
E* B — £"b cam to discriminate the signal from the back- 
ground, where £"b cam is the run-dependent beam energy, 
and Eg and p* B are the reconstructed energy and mo- 
mentum, respectively, of the B meson candidates in the 



center-of-mass (CM) frame. We retain B candidates 
with 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.15 
GeV < AE < 0.2 GeV. After all selection requirements, 
2.3% of the events for iP{2S){£+£-)ir- and 7.1% for 
iP(2S)(J/ijjtt + tt^)tt^ contain more than one B candidate. 
For these events, we choose the B candidate whose M\, c 
is closest to the nominal B meson mass. 

To suppress continuum background, events having a 
ratio of the second to zeroth Fox- Wolfram moments [21 1 
R 2 > 0.5 are rejected. Large B -> (J/ip,ip(2S))X MC 
samples corresponding to 3.86 x I0 10 generic BB events 
are used for the background study, because backgrounds 
predominantly come from B decays into final states hav- 
ing J/ip -> £+£- or ip{2S) -> £+£-. We find that 
the dominant background comes from B~ — > ip{2S)K~ 
where the kaon is misidentified as a pion. This decay 
mode makes a peak at AE ~ —0.07 GeV. Other back- 
grounds originate from random combinations of ip(2S) 
and n~ candidates (combinatorial background) and do 
not form any peaking structure in the AE projection of 
the MC sample. Studies of AM and £ + £~ invariant mass 
data sidebands support this assumption. 

We extract the signal yield by performing an unbinned 
extended maximum likelihood fit to the AE distribution 
of the selected B candidates. The extended likelihood 
function used is 



-(N S + N BK +N C ) N 

£(N s ,N BKl N c ) = — ^ Y[[N s P s (AE t y 

4=1 



Nl 



NsKPeKiAEi) + NcPciAE, 



(1) 



where N is the total number of candidate events; N$, 
Nbk and Nq denote the number of signal, background 
from B~ — ► ip(2S)K~ and combinatorial background 
events, respectively. The modeling PDF's include a sum 
of two Gaussians for the signal (Ps(AEi)), a sum of 
two bifurcated Gaussians for the B~ — > ip(2S)K~ back- 
ground (PBfc(AEi)) and a second-order polynomial for 
other backgrounds (Pc(AEi)). In this study the B~ — ► 
ip(2S)K~~ decay mode is used as a control sample, as 
well as the denominator of the relative branching frac- 
tion ratio B(B~ -> 4>{2S)n-)/B{B- -> ip(2S)K-). 
The B~ — > ip(2S)ir~ signal shape is fixed to that ob- 
tained from the B~ — >• ip(2S)K~ control sample. The 
B~ — > ijj(2S)K~ background shape is determined from 
the MC sample after applying a correction for the differ- 
ence between data and MC. For the smooth combinato- 
rial background, parameters of the second-order polyno- 
mial are floated in the fit. 

We perform a simultaneous fit to all the considered 
ip(2S) decay modes to obtain a common branching frac- 
tion, taking into account different detection efficiencies 
for each sub-decay mode and detector configuration. The 
efficiencies arc determined from signal MC samples after 
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FIG. 2: AE distributions of the B~ -> ip{2S)ir~ 
candidates, reconstructed in the following four sub- 
modes: (a) B~ -> il>{2S)(J/ip{ee)ir + n-)ir~, (b) B~ -» 
V'(2S')(J/^( W )7r+7r-)7r-, (c) B- -» V(25)(ee)7r- and (d) 
B~ — » ijj(2S)((j.(i)ir~ . The curves show the signal (green 
dashed) and the background components (red dot-dashed for 
B~ — > ip(2S)K~ , magenta dotted for combinatorial back- 
ground) as well as the overall fit (blue solid). 



applying the correction factors taking into account the 
data and MC differences as described later. In the fit, the 
yield of B~ — > ip{2S)K~ (background) and the yield and 
polynomial parameters of the combinatorial backgrounds 
are floated for each decay mode separately. The fit gives 
a branching fraction of (2.44 ±0.22 ±0.20) x 10~ 5 , where 
the first error is statistical and the second is systematic. 
Equal production of neutral and charged B meson pairs 
in the T(45) decay is assumed. We use the J/ip — > e + e~, 
J/ip -» n + n~, ip{2S) -> J/ipn+ir-, ip{2S) -» e+e- and 
ip{2S) — > fi + L l ~ intermediate branching fractions from 
Ref. [20j. Separate fits to the four sub-modes are also 
performed. The results are shown in Fig. 2 and sum- 
marized in Table I. The yields are calculated from the 
obtained branching fractions and efficiencies. The effi- 
ciencies listed in Table I are weighted averages of the two 
separate data sets. For all sub-modes, we observe signals 
with more than 5<r statistical significance and the branch- 
ing fractions agree well with one another. The statistical 
significance is defined as y/— 2 In {Co / C maK ) where C max 
{Cq) denotes the likelihood value at the maximum (with 
the signal yield fixed to zero). 

A correction for small differences in the signal de- 
tection efficiency calculated from signal MC and data 



has been applied for the pion identification requirement, 
muon identification and AM. Uncertainties on these cor- 
rections are included in the systematic error. The J/ip — > 
£ + £~ and e + e~ — > e + e~£ + £~ samples are used to esti- 
mate the lepton identification correction and uncertainty 
whereas the pion (kaon) identification correction and un- 
certainty are determined from a D* + — > D Q {K~ir + )ir + 
sample. A correction factor for the AM requirement is 
determined from the B~ — > ip{2S)K~ sample and is es- 
timated by taking the ratio of yields from data and MC 
for tight (0.578 GeV/c 2 < AM < 0.598 GcV/c 2 ) and 
loose (0.570 GcV/c 2 < AM < 0.610 GeV/c 2 ) selection 
cuts. The fitting procedure is checked using signal and 
B — ► {J/ip,ip{2S))X MC samples and no significant bias 
is found. The signal yield systematic uncertainty is calcu- 
lated by varying each fixed parameter in the fit by ±ler, 
and then taking the sum in quadrature of the deviations 
of the signal yield from the nominal value. The total sys- 
tematic uncertainty assigned to the yield estimation is 
3.5%. The uncertainties due to the daughter branching 
fractions amount to 3.4%. The uncertainty on the track- 
ing efficiency is estimated to be 1.2% per track and 5.0% 
in total. Systematic errors have been estimated by taking 
into account whether the particular error is correlated or 
uncorrelated and then combining them with the proper 
weights. All the systematic uncertainties are summarized 
in Table II. The total systematic error of 8.4% is the sum 
in quadrature of all the uncertainties. 

The CP-violating charge asymmetry A is defined as : 



A 



B{B- -> ip{2S)ir-) - B{B+ -> i/j{2S)it+) 
B{B- -* V(2S>-) + B{B+ -> V(25)tt+) ' 



(2) 



We extract branching fractions for B + and B~ samples 
separately using the same procedure as described above, 
except the polynomial background shape is fixed to that 
obtained in the previous fit to the total B^ — > V'(2S')7r ± 
yield. We obtain 



.4 = 0.022 ±0.085 ±0.016 



(3) 



with corresponding signal yields of 89 ± 11 and 93 ± 11 
for B + and B~ , respectively. Here the systematic error 
includes contributions from the uncertainty on the yield 
extraction method, charge asymmetry in the pion identi- 
fication efficiency as well as a possible detector bias and 
a possible difference between B + and B~ signal shape 
parameters. Among all these uncertainties, the only non- 
negligible one is a possible detector bias: it is estimated 
to be 0.016 from a study of B~ -> J/tpK~. 

We also measure the ratio of B{B~ — ► ip{2S)ir~) and 
B{B~ -> ip{2S)K~). In a study using large B -> {J/ip, 
ip{2S))X MC samples, the amount of background for 
B~ -> ijj{2S){£+£ )K- is found to be negligible. On the 
other hand, in the case of B~ — > ip{2S){J/ipir + ir~)K~ , 
the decay mode B~ —* J/ipKi{1270)~ as well as B~ — > 
J /ip'K Jr -K~ K~ appear as peaking backgrounds since they 



TABLE I: Summary of the results. 



Decay mode 


Efficiency(%) 


Signal yield 


Branching fraction(10 ) 


Statistical significance 


B~ -* iP(2S)(J/TP(e+e-)ir + iT-)iT- 
B~ -^ ^(2S)(J/^(fi + ^-)TT + n-)-K- 
B~ -^ip(2S)(e + e')n- 
B~ -> ^{2S)\ix + lT)-K~ 


15.1 
16.8 
32.2 

35.7 


48.9±8.3 
44.0±8.1 
44.0±9.0 
43.5±7.7 


2.57 ±0.44 
2.08 ±0.38 
2.80 ±0.57 
2.50 ± 0.44 


9.5cr 
8.4(7 
7.3cr 

9.0cr 


B~ -> tp(2S)TT~ (combined) 






2.44 ±0.22 ±0.20 





TABLE II: Summary of systematic errors on the B 
if)(2S)ir~ branching fraction. 



Source 



Uncertainty (%) 



Uncertainty in yield 
Tracking error 
Lepton identification 
Pion identification 
MC statistics 
Number of BB pairs 
Daughter branching fractions 
AM requirement 



3.5 
5.0 

4.2 
1.6 
0.3 
1.4 

3.4 
0.5 



Total 



8.4 



have the same final state as the signal. The background 
estimation is performed using the AM data sideband de- 
fined as 0.51 GeV/c 2 < AM< 0.57 GeV/c 2 and 0.61 
GeV/c 2 < AM< 0.67GeV/c 2 . We obtain 7.2±2. 7 events 
for B~ -> ip(2S)(J/tp(ee)nir)K- and 9.4 ±2.7 events for 
B~ —¥ ip(2S)(J/ip(fifi)irir)K~ when scaled to the sig- 
nal region. These backgrounds are subtracted from the 
appropriate signal yields. A simultaneous unbinned ex- 
tended maximum- likelihood fit gives a branching fraction 
of (6.12 ± 0.09) xlO -4 and corresponding signal yield of 
4720 ± 69, where the error is statistical only. The result- 
ing fit is shown in Fig. 3. The branching fraction is in 
good agreement with our previous measurement [4J and 
the world-average [20j. 

We obtain the ratio of B(B~ -> ip(2S)iv-) and 
B{B~ -* ij>(2S)K-) to be 



B(B~ -» i/)(2S)ir 



B(B- 



^(2S)K~) 



4- = (3.99 ±0.36 ±0.17)%, (4) 



which is consistent with the expectations of the factor- 
ization hypothesis. Many sources of systematic errors 
cancel in the ratio of the branching fractions. Contribu- 
tions to the systematic error come from the uncertainty in 
pion identification, signal extraction method, MC statis- 
tics from B~ — > ip(2S)Tr~ (K~) decay, kaon identifica- 
tion and uncertainty on the background estimation from 
B~ — > tp(2S)K~ decay. The total uncertainty is 4.2%. 

In summary, we report the first observation of the de- 
cay B~ -> ij)(2S)-K~~ using 657 x 10 6 BB events. The mea- 




FIG. 3: AE distribution of the B~ -> tp(2S)K~ candidates. 
The curves show the overall fit and the background compo- 
nent (magenta dotted for combinatorial background). 



sured branching fraction is B(B — > ip(2S)ir ) = (2.44 ± 
0.22 ± 0.20) x 10~ 5 . No significant direct CP- violating 
charge asymmetry is observed in B~ — ► t/j(2S)ir~ . The 
ratio of branching fractions is (3.99 ±0.36 ±0.17)% which 
is consistent with the theoretical prediction based on fac- 
torization 
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